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ABSTRACT 

The standard cold dark matter (CDM) model has recently been challenged by the 
claim that dwarf galaxies have dark matter haloes with constant density cores, whereas 
CDM predicts haloes with steeply cusped density distributions. Consequently, numer- 
ous alternative dark matter candidates have recently been proposed. In this paper, we 
scrutinize the observational evidence for the incongruity between dwarf galaxies and 
the CDM model. To this end, we analyze the rotation curves of 20 late-type dwarf 
galaxies studied by Swaters (1999). Taking the effects of beam-smearing and adiabatic 
contraction into account, we fit mass models to these rotation curves with dark matter 
haloes with different cusp slopes, ranging from constant density cores to r -2 cusps. 
Even though the effects of beam smearing are small for these data, the uncertainties 
in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density 
distribution hamper a unique mass decomposition. Consequently, the rotation curves 
in our sample cannot be used to discriminate between dark haloes with constant den- 
sity cores and r~ l cusps. We show that the dwarf galaxies analyzed here are consistent 
with cold dark matter haloes in a ACDM cosmology, and that there is thus no need to 
abandon the idea that dark matter is cold and collisionlcss. However, the data is also 
consistent with any alternative dark matter model that produces dark matter haloes 
with central cusps less steep than r -15 . In fact, we argue that based on existing HI 
rotation curves alone at best weak limits can be obtained on cosmological parameters 
and/or the nature of the dark matter. In order to make progress, rotation curves with 
higher spatial resolution and independent measurements of the mass-to-light ratio of 
the disk are required. 
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1 INTRODUCTION 

The standard cosmological model for structure formation 
combines an inflationary Universe with hierarchical growth 
of structures that originate from small fluctuations in the 
cosmic mass distribution. In addition to baryonic matter, 
this model requires non-baryonic dark matter and possibly 
some form of vacuum energy or quintessence. Unfortunately, 
the nature of the dark matter, which is the dominant mass 
component, still remains unknown. A large number of candi- 
dates have been proposed of which cold dark matter (CDM) 
has been the most popular. Because CDM particles have 
a negligible thermal velocity with respect to the Hubble 
flow, the original phase-space density of cold dark matter 
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is extremely high. Numerical simulations have shown that a 
small fraction of this material remembers its initial phase- 
space density even after it collapses to form a bound object. 
This low-entropy material settles in the centers of virialized 
dark haloes, thus creating steeply cusped density profiles, 
and causing a large fraction of haloes to survive as sub- 
structure inside larger haloes (e.g., Navarro, Frenk & White 
1996; Fukushige & Makino 1997; Moore et al. 1998, 1999b; 
Ghigna et al. 1998; Klypin et al. 1999a; White & Springel 
1999). 

These characteristics of CDM haloes, however, seem to 
disagree with a number of observations. First, the number 
of sub-haloes around a typical Milky Way galaxy, as iden- 
tified by satellite galaxies, is an order of magnitude smaller 
than predicted by CDM (Kauffmann, White & Guiderdoni 
1993; Klypin et al. 1999b; Moore et al. 1999a). Secondly, the 
observed rotation curves of dwarf and low surface bright- 
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ness (LSB) galaxies seem to indicate that their dark mat- 
ter haloes have constant density cores instead of steep cusps 
(Flores & Primack 1994; Moore 1994; Burkert 1995; Burkert 
& Silk 1997; McGaugh & de Blok 1998; Stil 1999; Moore et 
al. 1999b; Dalcanton & Bernstein 2000; Firmani et al. 2001). 
In view of these discrepancies, numerous alternatives to 
the CDM paradigm have recently been proposed. These in- 
clude broken scale-invariance (hereafter BSI; Kamionkowski 
& Liddle 2000; White & Croft 2000), warm dark matter 
(hereafter WDM; Sommer-Larsen & Dolgov 1999; Hogan & 
Dalcanton 2000), scalar field dark matter (hereafter SFDM; 
Peebles & Vilenkin 1999; Hu & Peebles 1999; Peebles 2000; 
Matos, Siddharta & Urena-Lopez 2000), and various sorts 
of self-interacting or annihilating dark matter (hereafter 
SIDM; Carlson, Machacek & Hall 1992; Spergel & Stein- 
hardt 2000; Mohapatra & Teplitz 2000; Firmani et al. 2000; 
Goodman 2000; Kaplinghat, Knox & Turner 2000; Bento 
et al. 2000). Whereas particle physics does not prefer CDM 
over for instance WDM, SFDM, or SIDM, the former has 
the advantage over the latter that it has no free parameters. 
Furthermore, most of these alternatives seem unable to solve 
both problems simultaneously (Moore et al. 1999b; Hogan 
& Dalcanton 2000; Coli'n, Avila-Reese & Valenzuela 2000; 
Dalcanton & Hogan 2000), and often the alternatives face 
their own problems (Spergel & Steinhardt 2000; Hannestad 
1999; Burkert 2000; Moore et al. 2000; Yoshida et al. 2000; 
Kochanek & White 2000; Miralde-Escude 2000; Sellwood 
2000). 

As an alternative to modifying the nature of the dark 
matter, the sub-structure and core problems might be solved 
once additional baryonic physics are taken into account. Sev- 
eral studies have suggested that processes such as reioniza- 
tion and supernova feedback can help to suppress star for- 
mation and to decrease central densities in low-mass dark 
matter haloes (e.g., Navarro, Eke & Frenk 1996; Gelato & 
Sommer-Larsen 1999; van den Bosch et al. 2000; Bullock, 
Kravtsov & Weinberg 2000; Binney, Gerhard & Silk 2001). 
Whereas these processes may indeed help to solve the prob- 
lem with the over-abundance of satellite galaxies, the sug- 
gestion that feedback processes can actually destroy steep 
central cusps seems somewhat contrived in light of more de- 
tailed simulations. For instance, as shown by Navarro, Eke 
& Frenk, the effects are only substantial if large fractions 
of baryonic mass are expelled, which seems hard to rec- 
oncile with the low ejection efficiencies found in more de- 
tailed hydro-dynamical simulations (e.g., Mac-Low & Fer- 
rara 1999; Strickland & Stevens 2000). 

It is evident from the above discussion that the long- 
time popular CDM paradigm is currently facing its biggest 
challenge to date. However, before abandoning CDM on the 
grounds that it is inconsistent with observations, it is worth- 
while to more closely examine the observational evidence 
against it. In this paper we scrutinize CDM's most persis- 
tent problem: the claim that the dark haloes of dwarf galax- 
ies, as inferred from their rotation curves, are inconsistent 
with CDM predictions. The main motivation for this work is 
that recent work on the rotation curves of LSB galaxies has 
shown that once data with sufficient resolution is obtained, 
or the effects of beam smearing are properly taken into ac- 
count, the inner rotation curves are significantly steeper and 
allow for more centrally concentrated dark matter haloes 
(Swaters 1999; van den Bosch et al. 2000; Swaters, Madore 



& Trewhella 2000). In fact, these studies have pointed out 
that, in contrast with previous claims, current data on LSB 
rotation curves are consistent with CDM predictions. 

Here we analyze a set of HI rotation curves of a sample 
of 20 late-type dwarf galaxies. Taking beam smearing and 
adiabatic contraction of the dark matter into account, we 
investigate whether the rotation curves of the galaxies in our 
sample are consistent with CDM haloes. Although we cannot 
rule out that these dwarf galaxies have dark haloes with 
constant density cores, we find that their rotation curves 
are consistent with cold dark matter haloes as expected in 
a ACDM cosmology. 



2 THE DATA 

The HI rotation curves that we use in this paper have been 
derived from the data presented in Swaters (1999, hereafter 
S99). The HI observations were done with the Westerbork 
Synthesis Radio Telescope (WSRT) as part of the West- 
erbork HI Survey of Spiral and Irregular Galaxies Project 
(WHISP, see Kamphuis, Sijbring & van Albada 1996). The 
observations and data reduction are discussed in detail in 
S99. From the sample of 73 late-type dwarf galaxies we only 
selected those galaxies that according to S99 have high qual- 
ity rotation curves (no strong asymmetries and high signal- 
to- noise ratio). This sample consists of 20 galaxies, which 
have inclination angles in the range 39° < i < 80° . Galaxies 
with i > 80° have been excluded from the sample because 
they require a somewhat different analysis. 

Table [I] lists the properties of the galaxies in our sam- 
ple. The absolute magnitudes, disk scale lengths, and cen- 
tral surface brightnesses have been determined from _R-band 
photometry (Swaters & Balcells 2001, hereafter SB01). The 
distances are as adopted by SB01: where possible stellar dis- 
tance indicators have been used, mostly brightest stars. If 
these were not available, distances based on group mem- 
bership were used, or, if these were unavailable as well, the 
distance was calculated from the HI systematic velocity fol- 
lowing the prescription in Kraan-Korteweg (1986) with an 
adopted Hubble constant of Ho = 75 kms -1 Mpc~\ B — R 
colors, available for only 6 galaxies, are also taken from 
SB01. 

The original HI observations have been obtained with a 
typical beam of 14" x 14"/ sin 5 (with 8 the object's declina- 
tion). In general, the signal-to-noise ratio at this resolution 
was too low to obtain reliable rotation curves. Therefore, 
the data were convolved to a resolution of approximately 
30" x 30". Velocity fields were constructed by fitting Gaus- 
sian curves to the observed line profiles at each position. 
Next, the rotation velocities and their formal errors were de- 
termined by fitting a tilted-ring model to the velocity fields 
assuming constant inclination and position angles. Where 
possible the orientation angles were determined from the 
velocity fields, and in the remaining cases from the optical 
images. For details about the determination of the tilted 
ring parameters see S99. Note, however, that our analysis 
is different from the one presented in S99, where an itera- 
tive method, based on modelling of the observed data cubes, 
was used to approximately correct the rotation curve for the 
effects of beam smearing. The rotation curves we use here, 
however, have not been corrected for beam smearing. In- 
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Table 1. Properties of sample of late- type dwarf galaxies. 
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Column (1) lists the UGC number of the galaxy. Columns (2) - (6) list the 
distance to the galaxy (in Mpc), absolute i?-band magnitude, central ij-band 
surface brightness (in mag arcsec - 2 ), scale length of the stellar disk (in kpc), 
and the observed rotation velocity Vi ast (in kms -1 ) at the last measured point. 
Column (7) lists the number of data points, Ny, along the rotation curve (two 
data points per beam). Columns (8) and (9) list the B — R color (if available) and 
the adopted inclination angle (in degrees), respectively. Magnitudes and central 
surface brightnesses have been corrected for inclination and galactic extinction, 
but not for internal extinction. Finally, column (10) indicates whether the galaxy 
is consistent with a ACDM cosmology (+) or whether no meaningful fit can be 
found (?). See the discussion in § H for details. 



stead, we beam-smear our models before comparison with 
the data, following the procedure detailed in § |3.2]. 



3 ROTATION CURVE FITTING 
3.1 Mass components 

For the mass modelling presented in this paper, we assume 
that there are three main mass components in each galaxy: 
an infinitesimally thin gas disk, a thick stellar disk, and a 
spherical dark halo. We closely follow the procedure outlined 
in van den Bosch et al. (2000; hereafter BRDB), which we 
briefly outline below for completeness. 

In order to determine the contribution of the gas to the 
galaxy's circular velocity, we make the assumption that the 
gas is distributed axisymmetrically in an infinitesimally thin 
disk. Under this assumption the circular velocities due to the 
self-gravity of the gas can be computed from the gas surface 
density using equation [2-146] of Binney & Tremaine (1987). 
We model the HI density distribution as follows: 



"(I) 



-R/Ri 



+ f So e 



-((fl-K 2 )/<r) 2 



(1) 



The first term is identical to the surface density profile used 
in BRDB, and represents an exponential disk with scale 
length Ri and with a central hole, the extent of which de- 
pends on /3. The second term corresponds to a Gaussian 
ring with radius R2 and a FWHM tx a. The flux ratio be- 
tween these two components is set by /. The form of equa- 



tion (|l|) has no particular physical motivation, but should be 
regarded as a fitting function. When computing the circular 
velocities of the atomic gas, we multiply Ehi by a factor 1.3 
to correct for the contribution of helium. 

For the stellar disk we assume a thick exponential 



p*(R,z) = po exp(— R/Rd) sech 2 (z/zo) 



(2) 



where Rd is the scale length of the disk. Throughout we set 
zo = Rd/6- The exact value of this ratio, however, does not 
significantly influence the results. The circular velocity of the 
stellar disk is computed using equation [A. 17] in Casertano 
(1983), and properly scaled with the stellar ii-band mass- 
to- light ratio Tr. None of the galaxies in our sample has a 
significant bulge component. 

We assume that initially the dark and baryonic matter 
virialize to form a spherical halo with a density distribution 
given by 



p(r) 



Po 



(r/r s ) a (l+r/r a ) 3 - 



(3) 



with r s being the scale radius of the halo, such that p tx r a 
for r <§; r s and p tx r~~ 3 for r ^> r s . For a — 1 equation ^ 
reduces to the NFW profile (Navarro, Frenk & White 1997). 
We define the concentration parameter c = r2oo/r s , with 
J"200 the radius inside of which the mean density is 200 times 
the critical density for closure, i.e., 



T200 V2OO 

/i _1 kpc kms -1 



(4) 
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Here V200 is the circular velocity at r2oo, and h = 
tfo/100 km s^Mpc -1 . 

The formation of the stellar and gaseous disks due to 
the cooling of the baryons inside the virialized halo leads to 
a contraction of the dark matter component. We make the 
assumption that the baryonic collapse is slow, and take this 
adiabatic contraction of the dark halo into account follow- 
ing the procedure in Barnes & White (1984), Blumenthal et 
al. (1986) and Flores et al. (1993). The halo mass inside ra- 
dius r, required for the adiabatic contraction computations, 
is given by 



M halo (r) = M200 



fi(xc) 



with x — r/r 200 and 



H(x) 



- a (l + y) a - 3 dy 



(5) 



(6) 



3.2 Beam smearing 

As is evident from the results presented in S99 and BRDB, it 
is important that the effects of beam-smearing are properly 
taken into account. Rather than attempting to deconvolve 
the observations (which is an ill-defined problem), we con- 
volve our models with the effective point spread function P 
(i.e., the beam) of the interferometer. The convolved surface 
brightness at a position (x, y) on the plane of the sky is 



T,{x,y) 



drr dOE(r')P(r,0-O o ). 



(7) 



Here r' — -J x' 2 + y' 2 , where x = x + rcos9 and y' = 
(y + r sin 9) / cos i are the Cartesian coordinates in the equa- 
torial plane of the disk, i is the disk's inclination angle, and 
9o is the angle between the major axes of the galaxy and the 
beam (for which we adopt a two-dimensional Gaussian, see 
BRDB). The underlying surface brightness, E(7?), is mod- 
eled by equation (|l|). Note that we assume that the HI dis- 
tribution is axisymmetric. 

Beam smearing also affects the observed rotation ve- 
locities Vrot at a position (x,y) on the plane of the sky by 
causing gas from a larger area of the disk to contribute to 
the observed line of sight velocity: 



V ro t(x,y) = ^ / drr d9 E(r')Vi oa (s;', y')P(r, 9 — 6>o) (8) 



where Vi os is the line of sight velocity. Throughout we as- 
sume that the gas moves on circular orbits in the plane of 
the disk and has zero intrinsic velocity dispersion. As ar- 
gued by Swaters (1999), the asymmetric drift corrections as 
calculated from the HI distribution are generally small and 
have little effect on the derived rotation curves. 



3.3 Fitting procedure 

The first step in fitting mass models to the rotation curves 
is to determine the best-fit model for the true underlying 



HI distribution. The surface density distribution of equa- 
tion (|l|) has six free parameters. Note, however, that So is 
completely determined by normalizing the models to the to- 
tal mass in HI and can thus be ignored in the fitting routine. 
We determine the best-fit parameters by minimizing 



2 

Xhi = 



JVhi 

E 



(9) 



with E bs the observed HI density distribution at 30" reso- 
lution and AE bs the corresponding errors. The resu lts are 
shown in the upper- right panels of Figures A1-A19. Open 



circles correspond to the observed HI surface density and 
solid lines to the best-fit model. In most cases, the un- 
smeared HI distribution of the best-fit models, indicated 
by dashed lines, is fairly similar to that after convolution 
with the beam, indicating that the effects of beam-smearing 
for these data are only small. In two cases, UGC 7524 and 
UGC 7603, our fitting function (equation fl|) can not satis- 
factorily describe the data. In these cases we use the data of 
the full resolution (see § []) as a model for the true underlying 
HI distribution. 

Once Ehi is known we can compute the beam-smeared 
model rotation curves. For a given choice of the Hubble con- 
stant the mass models described above have four free param- 
eters to fit the data: Tr, a, c, and V200 (or equivalently r2oo). 
For a given («,Th) we determine the best-fitting c and V200 
by minimizing 



2 

Xvei 



i=l 



AVobs(ft) 



(10) 



Here AVobs are the formal errors on Vobs from the tilted-ring 
model fits, and V(Ri) is computed from equation ([j) with 
x = Ri and y — 0. 



3.4 Uncertainties on the rotation velocities 

An important issue in constraining the density distribution 
of the dark matter haloes is how to interpret Xvei- O ne can 
only use the absolute values of Xvei 1° compute confidence 
levels for our models, if the errors AVobs are the proper, 
normally distributed errors, there are no systematic errors, 
the data points are independent, the assumptions underlying 
the model are correct, and the mass-model is a proper rep- 
resentation of the real mass distribution. However, the fact 
that the minimum x V ci differs considerably from the num- 
ber of degrees of freedom for almost all galaxies indicates 
that we do not meet these criteria. This does not come as a 
surprise. First of all, errors in the assumed inclination angle, 
distance, beam parameters, and the distribution of gas and 
stars all lead to systematic errors in the dark matter density 
distribution. Furthermore, there are numerous assumptions 
underlying our mass-models, each of which may be in error. 
For instance, we assume that the halo is sperical, that the 
disk is axisymmetric and that the gas moves on perfectly 
circular orbits with zero intrinsic velocity dispersion (i.e., 
we thus ignore asymmetric drift). In addition, we assume 
that E gas = 1.3Ehi and that T_r is constant throughout the 
stellar disk. We thus ignore any contribution from molecular 
and/or ionized gas as well as any radial changes in stellar 
population. Given all these potential sources of confusion, 
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we only use Xvei to assess the relative quality of the model 
fits. We do not try to assign any confidence levels to the 
absolute values of either \vei or ^Xvei- Henceforth, if, for in- 
stance, a model for a particular galaxy yields a smaller Xvei 
for a — than for a = 1, the model with a constant density 
core provides a better fit than the NFW model, but it does 
not necessarily imply that the rotation curve is inconsistent 
with CDM. For these reasons, we use physical criteria (dis- 
cussed in § s|) rather than criteria based on uncertain con- 
fidence levels to assess whether a model is consistent with 
the data. 



3.5 Degeneracies in the mass modelling 

Before interpreting the results in terms of constraints on 
the density distribution of dark matter haloes it is useful 
to examine some models. To that end, we construct three 
model galaxies moulded after UGC 731, i.e., the models have 
the same HI and stellar disks as UGC 731, and we adopt 
the same distance and inclination angle (see Table |^). We 
add a dark matter halo with a density distribution given by 
equation (E3) and compute the beam-smeared rotation veloc- 
ities at 15 intervals. This corresponds to roughly half the 
beam size and is identical to the interval between actual data 
points of the rotation curve of UGC 731. We convolve the 
velocity field using the same beam-size and beam orienta- 
tion as for the true UGC 731 data. Finally, we add a random 
Gaussian error to the model rotation velocities (with vari- 
ance AV). All three models have the same mass distribution: 
a = 1, c = 20, V200 = 75 kms -1 , and T R = 2.0 (M/L) . 
What we vary, however, is the way the rotation curve is sam- 
pled: Model 1 has a rotation curve with 6 data points ex- 
tending out to 0.06 r2oo and with AV = 2.0 kms -1 . The ro- 
tation curve of model 2 has twice as many data points, thus 
extending twice as far out, and has the same AV. Model 3, 
finally, has a rotation curve that extends equally far as that 
of model 2, but with AV — 0.2 kms -1 . For comparison, 
the real data of UGC 731 consists of 12 data points with 
(AV) — 1.4 kms , and model 2 is thus a fair representa- 
tion of the actual data. Also, as can be seen in column (7) 
in Table [lj the number of data points for the other galaxies 
in our sample varies from 5 (UGC 7232) to 31 (UGC 7524). 
Note that, because of distance effects, a larger Nv does not 
necessarily mean that the density distribution of the halo is 
probed to larger radii. 

We analyze our model rotation curves in the same way 
as we analyze the data for the dwarf galaxies. The results 
are shown in Figure [j]. A number of general tre nds, w hich 
are also present in the real data (see Figures A.1-A1E), are 
immediately apparent. First of all, the halo concentration 
c of the best-fit model decreases with increasing Tr and 
a. This is easily understood in terms of the total enclosed 
mass which has to be similar for different models. Secondly, 
there is an a cr i t for which V200 is maximal and c = 1. For 
a > Qcrit the best fitting halo concentration c < 1, implying 
that the scale radius r s is larger than the virial radius r2oo- 
Clearly, for c < 1, equation (H) no longer is an appropriate 
description of the dark matter distribution, and we there- 
fore demand that c > 1. Consequently, for a > a cr it the 
halo concentration c = 1 and the quality of the fits rapildy 
decreases (i.e., \ 2 increases). 

In addition to these trends, a large amount of non- 



uniqueness is apparent, which can be associated with two 
distinct degeneracies. First of all, for any given cusp slope a, 
the relative amount of mass in the stellar disk can be traded 
off against the amount of mass in the halo, while maintaining 
virtually equally good fits to the data, i.e., different com- 
binations of (Tr,c, V200) yield similar values of xlei- This 
degeneracy is well-known from the rotation curves of high 
surface brightness galaxies (e.g., van Albada et al. 1985), and 
is generally referred to as the mass-to-light ratio degeneracy. 
The second degeneracy is that for a given mass-to-light ra- 
tio Tn, different combinations of (a,c, V200) with a 1.5 
yield virtually equally good fits to the data, unless the errors 
on the observed rotation velocities are sufficiently small (cf. 
models 2 and 3). This degeneracy, which we refer to as the 
cusp-core degeneracy, is a consequence of the fact that rota- 
tion curves only sample a small fraction of the dark matter 
density distribution (see also Lake & Feinswog 1989). 

To illustrate the nature of this cusp-core degeneracy 
we construct the circular velocity curve, V a =i(r), of a dark 
matter halo with a = 1, c = 25, and V200 = 100 kms -1 . 
Next, for a range of values for a, we seek the values of c and 
V200 for which V a (r) (with a 7^ 1.0) best fits V a =i(r) out 
to a certain radius r max . The results are shown in Figure [| 
for r max = 0.15 r2oo (indicated by a dotted vertical line). 
The thick curve in the upper panels corresponds to V a =i(r), 
normalized to Vboo- The thin curves correspond to the best- 
fitting V a (r) for a = 0.0,0.2,0.4, 1.8. As is evident from 
the left panels in Figure ^, where we plot the circular ve- 
locities only out to 0.2r2oo, the different V a (r) curves are 
in fact very similar. Only for a <; 1.5 does V a (r) start to 
deviate more significantly from V a =i(r). This explains why 
the reduced Xve\ °f models 1 and 2 increases strongly for 
a <; 1.5. For a ^ 1.5 the circular velocity curves out to 
r max are remarkably similar, and only very accurate rota- 
tion curves (i.e., with small AV) can discriminate between 
the different curves (cf. models 2 and 3). Alternatively, ac- 
curate constraints on the actual density distribution of the 
dark matter requires a rotation curve that either extends 
sufficiently far, or that has sufficient independent measure- 
ments at very small radii. This is evident from the lower 
two panels of Figure ^| which plot the normalized difference 
(V a — U a= i)/V200 as function of the normalized radius, and 
which show how the different V a curves diverge for both 
r/r2oo <; 0.2 and r/r2oo 0.02. Unfortunately, in practice 
HI rotation curves rarely extend to large enough radii, do 
not have enough spatial resolution, have too large errors, and 
suffer too much from systematic effects to lift this cusp-core 
degeneracy. 



4 RESULTS 

The results for each individual galaxy are presented in Fig- 



ures Al 



Alt and discussed in some detail in the Appendix. 



For UGC 7557 no converging fit to the observed rotation 
curve could be achieved with the mass-model described in 
Section ^, and therefore no results are plotted for this galaxy 
(see discussio n in the Appendix). The four left panels of Fig- 



Al 



A19 show (from top to bottom) \r 

(see 



, c, V200, and 

the corresponding baryon fraction /bar (see § 3.5), all as a 
function of the cusp slope a, and for three different values 



of T R . Here Xr. 



Xvel/Ndl, with Ndi the number of de- 
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MODEL 1 MODEL 2 MODEL 3 




0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 



a a a 

Figure 1. Results of analyzing three model rotation curves, moulded after UGC 731. All models have the same density distribution 
(indicated by a black dot); they only differ in the extent of the rotation curve and the errors on the rotation velocities, as indicated in 
the text. In fitting the model rotation curves three different mass-to-light ratios have been assumed: T^j = 0, Tr = 2.0 (M/L)q (which 
is the input value of the models), and Tr = 4.0 (M/L)©. The labeling is as indicated in the left middle panel. Note that the mass 
models derived from the rotation curves of models 1 and 2 (both with AV = 2 kms' 1 ) are strongly degenerate. The rotation curve of 
model 3, for which AV = 0.2 kms -1 , however, allows a fairly accurate recovery of the input mass model, although the mass-to-light 
ratio degeneracy remains. 



grees of freedom. The same general trends we found for the 
model galaxies are also apparent in the data, i.e., the halo 
concentration c decreases with increasing Tr and a, and 
above a certain value of a the quality of the fits decreases 
rapidly while c = 1. As we indicated in § ^ we have used data 
that is convolved to a lower resolution in order to enhance 
the signal-to-noise ratio. We have checked that our best fit 
models are consistent with the higher resolution data and 
found good agreement. 



3.4 



above, we can not sim- 



As we have argued in 
ply use Xicd to put confidence levels on the various mod- 
els. Furthermore, we have pointed out that two distinct de- 
generacies hamper a unique mass-decomposition, which is 
readily apparent from the fact that models with very dif- 
ferent cusp slopes and/or mass-to- light ratios yield roughly 
equally good fits (see for instance UGC 11707, UGC 12060, 
and UGC 12632). However, some constraints can be im- 
posed by only considering models that are physically re- 



alistic. For instance, models with Tr = are unrealistic, 
and are therefore not considered meaningful model fits. In 
addition, we can use the baryon fraction of each model to 
check its physical validity. For each model we compute the 
baryonic mass fraction /bar = {M gas + M a tars)/Af200 with 
M200 = r2ooV£ 00 /G the total mass of the galaxy (baryons 
plus dark matter). For currently popular cosmologies with 
flo = 0.3 and h — 0.7, and using recent Big Bang nucleosyn- 
thesis constraints on the baryon density (fibar = 0.02/i -2 ; 
Buries & Tytler 1998), one expects a universal baryons frac- 
tion of roughly 0.14. Note that because we ignore any molec- 
ular and/or ionized gas, and because feedback may drive 
galactic winds and expel baryons from the halo, /bar may 
be significantly lower than the universal value. In order to 
leave some room for the uncertainties in the cosmological 
parameters, in what follows we shall only consider models 
unrealistic if /bar > 0.2. In section 5 below we address the 
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Table 2. Best-fit parameters for models with a = 1. 



0.4 0.6 

r / r a» 



Figure 2. An illustration of the cusp-core degeneracy. The thick 
solid lines in the upper two panels correspond to the circular ve- 
locity of a halo with a density distribution of equation (^) with 
a = 1.0 and c = 25. The abscissa and ordinate are normalized 
to r200 and V200 of this density distribution, respectively, The 
thin lines are best-fits to the inner parts of this circular veloc- 
ity curve of models with a = 0.0, 0.2, 0.4, 1.8, once again nor- 
malized to r2oo and V200 of the a = 1.0 model. For clarity, the 
model with a = 0.0 is indicated by dashed lines. The lower pan- 
els plot the normalized differences, (V a — V a =\) /V20Q as function 
of r/r200- When fitting the models, only the velocities out to 
r m „ = 0.15r2oo (indicated by vertical dotted lines) are taken 
into account. 



effects a possible lower bound on /bat might have on our 
conclusions. 

We can put some further constraints on the mass mod- 
els by considering what range of i?-band mass-to-light ra- 
tios to expect. Using the Bruzual & Chariot (1993) stellar 
population models we have computed B — R and Tr for a 
Scalo (1986) IMF and two different star formation histories. 
Figure ^ plots Tr versus B — R for three different metallici- 
ties and for both a constant star formation rate (left panel) 
and a single burst stellar population (right panel). For the 
six galaxies in our sample for which SB01 obtained B-band 
photometry we find {B - R) = 0.87 ± 0.09 (see Table [!]). 
If we assume that there is no internal extinction in these 
galaxies, this implies 0.5 5s Tr 1.1 for the stellar popula- 
tion models investigated here. This ignores the contribution 
of any non-luminous baryonic component that may have the 
same radial distribution as the stars, and which would in- 
crease the mass-to-light ratio. In Table ^ we list the param- 
eters of the best-fit models with a = 1 for both T_r = 
and Tr = 1.0 (M/L)q. Although models with Tr = are 
unrealistic, these best-fit models yield a useful upper limit 
on c (cf. Pickering et al. 1997; Navarro 1998). The mod- 
els with T_r = 1.0 (M/L)0 are chosen to represent a typ- 
ical mass-to-light ratio. Furthermore, a comparison of two 
models with different mass-to-light ratios sheds light on the 
(non)-uniqueness of the mass models. 

The panels in the lower-right corners of Figures |Al| - 



UGC 
731 

3371 

4325 

4499 

5414 

6446 

7232 

7323 

7399 

7524 

7559 

7603 

8490 

9211 

11707 

11861 

12060 

12632 

12732 



0.0 

1.0 

0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 
0.0 
1.0 



16.0 
13.5 
9.5 
8.0 
30.9 
25.7 
9.0 
1.6 

< 1.0 

< 1.0 
17.4 

7.0 
4.2 

< 1.0 
4.7 

< 1.0 
23.1 
15.2 

8.5 
4.8 
1.4 
1.2 
5.5 

< 1.0 
24.2 
13.5 
18.3 
14.8 
13.1 
11.2 
16.0 
12.5 
34.6 
24.3 
14.0 
12.6 
10.3 

6.8 



V200 
51.3 
52.3 
68.6 
69.8 
53.5 
52.7 
58.1 
131.5 
253.6 
128.8 
52.0 
59.5 
116.0 
59.6 
129.3 
193.9 
62.8 
70.9 
71.9 
82.5 
135.4 
135.1 
87.9 
217.3 
50.3 
57.2 
43.6 
44.5 
67.3 
66.9 
106.3 
100.6 
46.6 
46.5 
51.6 
51.4 
67.6 
73.3 



4 
3 
2 
2 
-4 



3 
-2 



/"bar 

2.3 x 10~ 2 
2.8 x 10~ 2 

1.6 x lO' 2 
2.1 x 10~ 2 
2.1 x 10~ 2 
4.1 x 10- 2 

2.7 x 10~ 2 

3.7 x 10" 3 

1.8 x 10 
2.3 x 10 
4.1 x 10 
4.8 x 10 

2.3 x 10 

3.8 x 10~ 3 

1.4 x 10~ 3 

1.5 x 10 

1.3 x 10 

1.4 x 10^2 
1.1 x 10~ 2 
1.4 x 10^2 
1.4 x 10~ 4 

1.6 x 10~ 4 

3.0 x 10~ 3 

3.9 x 10- 4 
3.3 x 10^2 
2.9 x 10^2 

5.6 x 10^2 

6.1 x 10~2 
5.3 x 10-2 
7.1 x 10-2 

x 10-2 
x 10-2 
7.8 x 10-2 
1.1 x 10- 1 
2.8 x 10-2 

3.7 x 10-2 
5.3 x 10-2 
4.7 x 10-2 



2.0 
7.0 



4.7 
4.7 
4.4 
4.4 
4.6 
4.6 
4.6 
3.7 
3.0 
3.7 
4.7 
4.5 
3.8 
4.5 
3.7 
3.3 
4.5 
4.4 
4.3 
4.2 
3.7 
3.7 
4.1 
3.2 
4.7 
4.6 
4.9 
4.8 
1.4 
1.4 
3.9 
4.0 
4.8 
4.8 
4.7 
4.7 
4.4 
4.3 



(c) 
15.6 
15.6 
14.6 
14.5 
15.5 
15.5 
15.2 
12.3 
10.1 
12.4 
15.6 
15.1 
12.8 
15.1 
12.4 
11.0 
14.9 
14.5 
14.4 
13.9 
12.2 
12.2 
13.7 
10.6 
15.6 
15.2 
16.2 
16.1 
14.7 
14.7 
13.1 
13.2 
16.0 
16.0 
15.6 
15.6 
14.7 
14.3 



31.3 
31.1 
29.2 
29.1 
30.9 
30.9 
30.4 
24.7 
20.2 
24.9 
31.1 
30.1 
25.5 
30.1 
24.9 
22.0 
29.8 
29.0 
28.9 
27.9 
24.4 
24.4 
27.4 
21.2 
31.3 
30.5 
32.4 
32.2 
29.3 
29.3 
26.1 
26.4 
31.9 
31.9 
31.1 
31.1 
29.3 
28.7 



For each galaxy we list the parameters of two best-fit models with 
a = 1 (i.e., a dark matter halo with a r _1 density cusp): one 
with Yr = and the other with = 1.0 (M/L)q. In addition 
to the best-fit parameters c, V200 (in kms -1 ), and the resulting 
baryon fraction /b ar , we list the mean halo concentration, (c) , for 
ACDM haloes with the same V200 as the best-fit model, as well 
as the 2<r lower and upper limits (c m ; n and c max , respectively) 
of the distribution of c. These values are computed using the 
model of Bullock et al. (1999) for the ACDM cosmology used 
here. Consistency with ACDM requires that c m j n < c < c max 
and / bar < 0.2. 



lISL plot the observed rotation curves (open circles with 
errorbars) together with four models with Tr — l.O(M/L)0. 
These are the best-fit models for a = (solid lines), a = 
0.5 (dotted lines), a = 1.0 (short-dash lines), and a — 1.5 
(long-dash lines). These plots illustrate the typical quality 
of the model fits and the dependence on the cusp slope a. 
In most cases the individual curves for the four models can 
hardly been discerned, further emphasizing the cusp-core 
degeneracy discussed above. 
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Figure 3. The i?-band mass-to-light ratio T^j as function of the B — R color of Bruzual & Chariot (1993) stellar population models with 
a Scalo IMF and either a constant star formation rate (left panel) or a single burst of star formation (right panel). Results are shown for 
three different metallicities: one fifth Solar (dotted lines), Solar (solid lines) and two and half times Solar (dashed lines). For the dwarf 
galaxies in our sample with known B — R we find (B — R) = 0.87 ± 0.09, which implies that ~ 0.5 — 1.1 (M/L)©. 



5 IMPLICATIONS FOR THE NATURE OF 
THE DARK MATTER 

The main goal of this paper is to assess whether or not 
the rotation curves of dwarf galaxies are consistent with 
CDM. Before we can address this, we need to define what 
"consistent with CDM" means in the context of the den- 
sity distribution of dark haloes. Ideally, one would like to 
place confidence levels on whether or not a rotation curve 
is consistent with certain CDM predictions. However, as we 
discussed in § i.4, we can not use our \ 2 statistic to com- 
pute such confidence levels. Therefore we follow a different 
approach. Several studies that have claimed inconsistencies 
between (dwarf) galaxy rotation curves and CDM haloes, 
indicated that when a model with a = 1.0 is fitted the im- 
plied halo concentration c is too low to be consistent with 
CDM (i.e., Pickering et al. 1997; Navarro 1998). We follow 
this approach and investigate whether the distribution of c 
for our best-fitting models with a = 1.0 is consistent with 
predictions for a particular CDM model. 

High resolution iV-body simulations have shown that 
haloes virialize to density distributions of the form of equa- 
tion (|J) with a ~ 1. Different simulations, however, often 
yield different values for the concentrations. Furthermore, 
the distribution of halo concentrations is fairly broad, and 
its median depends on the mass of the halo, its redshift, and 
the cosmology (Cole & Lacey 1996; Navarro, Frank & White 
1996, 1997; Avila-Reese et al. 1999; Bullock et al. 1999; Jing 
2000; Jing & Suto 2000). Henceforth, there is (currently) 
no well-defined boundary for c to be considered "consistent 
with CDM" . Instead, one can only ask whether the statisti- 
cal properties of a sample of rotation curves are consistent 
with CDM for a given cosmology and according to a given 
set of simulations. 

In what follows we focus on the currently popular 



ACDM cosmology with Q = 0.3, Qa = 0.7, h = 0.7, and 
as = 1.0. The particular simulations to which we compare 
our results are presented in Bullock et al. (1999; hereafter 
B99), who also give a simple recipe for computing c and its 
scatter as function of halo massQ. The reason for using this 
particular recipe is that it is tested against a large statisti- 
cal sample of several thousand haloes. The expected value 
of c for each galaxy (as determined by V200 of the best-fit 
model), and the 2a deviations from the median, indicated 
by c m i n and c max , are listed in Table |i[ 

We now define a rotation curve to be consistent with 
ACDM if for a = 1.0 and 0.5 <T R < 1.1 there is a best-fit 
model with c m i n < c < c max and /bar < 0.2. These limits 
on Tr and /b ar are motivated in § ^j. Note that we do not 
demand that the models with a = 1.0 yield the best fit (i.e., 
minimum Xred) °f all models. Instead, we demand that there 
is a model with a = 1.0 for which the resulting parameters 
are realistic. According to this definition, we find that 14 out 
of the 20 dwarf galaxies in our sample are consistent with 
ACDM. These galaxies are indicated by a '+' in column (10) 
in Table |l[ 

But what about the other 6 galaxies? For UGC 7577 
no converging fit could be obtained at all, and this galaxy is 
discarded from the following discussion. The best fit models 
of the other five galaxies all have a — 0.0 (constant den- 
sity core) and T_r = (which is unphysical). In Figure [| 
we plot /bar and V200 for these five galaxies as functions of 
Tr, whereby a is kept constant at zero. As can be seen, the 
inferred values of /bar rapidly decrease with increasing Tr, 
while V200 rapidly increases. For T_r <^ 0.5 Mq/ Lq, which 



t This model uses somewhat different definitions for the halo 
mass and concentration, which we convert to the c and M200 
used in our analysis. 
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Figure 4. The baryon fractions (left) and virial velocities (right) as function of Yr for the five galaxies that are inconsistent with ACDM 
haloes: UGC 5414 (solid lines), UGC 7232 (dotted lines), UGC 7323 (short-dashed lines), UGC 7559 (long-dashed lines) and UGC 7603 
(dot-dashed lines). The results plotted here are for a = (haloes with constant density cores) only. Note how /b ar deceases rapidly 
with increasing T^, while V200 rapidly increases. For T^ > 0.5 Mq/ Lq (which is the more realistic regime), /b ar < 8 X 10 -3 , while 
V200 > 90kms — 1 (except for UGC 7559). Given our current understanding of feedback, such low baryon fractions are hard to reconcile 
with the large halo masses, and we therefore consider these rotation curve fits unrealistic (see discussion in text). 



is the more realistic regime, one finds values of /bar well be- 
low 0.01. With a universal baryon fraction of ~ 0.14 this 
implies that > 90 percent of the available baryons would 
have to be expelled from the disk. Although current under- 
standing of feedback is very limited, and we are reluctant 
to impose stringent lower limits on /bar, such high ejection 
efficiencies seem inconsistent with the high values of V200 in- 
ferred from the best-fit models (see e.g., Dekel & Silk 1986; 
Efstathiou 2000). In particular, recent hydro-dynamical sim- 
ulations have indicated that even starburst driven winds 
are extremely inefficient in expelling matter from systems 
with such high virial velocities (Mac Low & Ferrara 1999; 
Strickland & Stevens 2000). Thus, unless a mechanism can 
be devised that can expell <; 90 percent of the available 
baryons from haloes with V200 ^ 90 kms -1 , we conclude 
that for these five galaxies none of the (a,T_R)-models are 
realistic. This implies that whereas the rotation curves of 
these galaxies are inconsistent with the ACDM model, they 
do not support an alternative picture in which, for instance, 
dark matter haloes have constant density cores. For these 
galaxies either (1) our mass-model is inadequate, (2) there 
are systematic errors i n the data, or (3) one or more of the 



assumptions listed in § 3.5 are wrong. 

The good agreement between CDM predictions and the 
majority of the rotation curves analyzed here is also illus- 
trated in Figure ^ where we plot c as function of V200 for 
the best-fit models with a = (left panels), a = 1.0 (mid- 
dle panels; see also Table ^), and a = 1.5 (right panels). 
Results are plotted for both Tr — (upper panels) and 
Tjj = 1.0 (M/L)q (lower panels). Solid circles correspond 
to galaxies that are consistent with ACDM, open circles 
to galaxies for which no meaningful fit can be obtained. 
The solid and dashed lines in the middle panels indicate 



the mean and the 2a intervals of the predictions based on 
the B99 model. For comparison, we also plot (dotted lines) 
the predictions based on the model of Navarro, Frenk & 
White (1997)0. For a = 0.0 and a = 1.5 no such lines are 
plotted, since no model predictions exist for these cases. As 
already noted by Eke, Navarro & Steinmetz (2000), the B99 
model predicts concentration values a factor ~ 1.6 larger at 
V200 = 30 kms -1 , and a slightly steeper mass dependence. 
Since the disagreement between the two models is relatively 
small compared to the expected scatter in c, our conclusions 
are not sensitive to the fact that we focus our discussion on 
the B99 model. It is apparent that for all galaxies for which 
a meaningful fit is obtained, the best fit values for c and 
V200 are consistent with the expected values in a ACDM 
cosmology. 

In total we thus find that 14 out of 20 galaxies are con- 
sistent with CDM. For the remaining 6 galaxies no mean- 
ingful fit to the observed rotation curves can be obtained 
with our mass models for any value of a. Henceforth, these 
galaxies are neither consistent with CDM, nor with any 
other viable alternative. We thus conclude that at present 
there is no convincing evidence that dwarf galaxies (or low 
surface brightness galaxies, see BRDB) have dark matter 
haloes that are inconsistent with CDM. However, we wish 
to point out that this conclusion is based on the presump- 
tion that CDM haloes have a ~ 1.0. If future simulations 
confirm the results by Fukushige & Makino (1997), Moore 
et al. (1998) and Klypin et al. (2000), that CDM produces 
more steeply cusped dark matter haloes with a ~ 1.5, we 



t Computed using the procedure outlined in their Appendix with 
/ = 0.01 and C = 3.41 X 10 3 
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Figure 5. The logarithm of the best-fit concentration parameter as function of the logarithm of V200 of the best-fit model. Results are 
shown for two different stellar mass-to-light ratios and three different values of a (as indicated in the panels). Solid circles correspond to 
galaxies that are consistent with the ACDM model, whereas open circles indicate galaxies for which no meaningful fit can be obtained 
(indicated by '+' and '?' in Table [l| respectively). The solid and dashed lines in the middle panels indicate the mean and the la limits 
of the distribution of halo concentrations as predicted by the B99 model for the ACDM cosmology. All galaxies for which the best-fit 
models are physically realistic are consistent with this model. However, for a = 1 .5 (right panels) it is apparent that a significant fraction 
of the best-fit models are unrealistic in that they have c = 1). Henceforth, if future simulations confirm that CDM yields haloes with 
a ~ 1.5 rather than a ~ 1.0, the rotation curves analyzed here may signal a true problem for the CDM paradigm. 



would have to conclude that the rotation curves analyzed 
here are only marginally consisten t wit h CDM. This is evi- 
dent from the plots in Figures A1-A19 and Figure 0, which 



show that for a <; 1.5 the fits in general become rather 
poor, and often unrealistic. Furthermore, we have made 
the assumption, based on stellar population models, that 
T_r 1.1 (M/L)©. If, however, the true mass-to- light ratio 
of the stellar disks is significantly higher and/or dwarf galax- 
ies have large amounts of (centrally concentrated) molecular 
gas, some of the galaxies analyzed here will have dark matter 
haloes that are inconsistent with CDM. 

Although we have shown that the majority of the rota- 
tion curves in our sample are consistent with CDM, it does 
not imply that they are inconsistent with any of CDM's al- 
ternatives, such as WDM, SFDM or SIDM. In fact, in some 
cases models with a constant density core provide better fits 
to the rotation curves than for a = 1.0 (e.g., UGC 7524 and 
UGC 9211). Ultimately, one might hope to use the rota- 
tion curves of (dwarf) galaxies to put some constraints on 
the central phase space densities or core radii of their dark 
matter haloes. This in turn constrains the masses and/or 



interaction cross sections of the dark matter particles. Un- 
fortunately, the rotation curves analyzed here do not put 
any significant constraints on the actual nature of the dark 
matter: they are consistent with any dark matter species 
that yields haloes with ^ a ^ 1.5. In order for the rota- 
tion curves to put stringent constraints on the nature of the 
dark matter, we have to be able to much better constrain 
the density distribution of the dark matter haloes. However, 
given the numerous potential sources of systematic errors 
and both the mass-to-light ratio degeneracy and the cusp- 
core degeneracy discussed above it seems unlikely that HI 
rotation curves alone will be able to provide any significant 
constraints on the nature of the dark matter. 



5.1 Comparison with previous work 

Our conclusion that the rotation curves of dwarf galaxies are 
consistent with a ACDM cosmology is at odds with previous 
studies (Moore 1994; Navarro et al. 1996; Moore et al. 1999b; 
Blais-Ouellette, Amram & Carignan 2000) . There are several 
reasons that may contribute to this discrepancy. First of all, 
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the previous studies have mostly used the same, small sam- 
ple of dwarf galaxies, the rotation curves of which have been 
determined in different ways. We have used high-resolution 
HI rotation curves for a new, relatively large sample of dwarf 
galaxies for which the rotation curves have been derived in a 
uniform way. Furtermore, some of the dwarf galaxies in the 
previous studies have high inclinations or strong non-circular 
motions e.g., as the result of a bar. Rotation curves derived 
from such galaxies are likely to suffer from systematic ef- 
fects. For this reason, we excluded galaxies with strong bars 
and inclinations larger than 80°. In addition, we have im- 
proved upon previous studies by taking beam-smearing and 
adiabatic contraction into account. As beam-smearing can 
mimic the presence of a constant density core, it is imper- 
ative that these effects are properly accounted for (see e.g., 
BRDB). Also, we have stressed the importance of degenera- 
cies, systematic errors, and the fact that not all models are 
physically meaningful. Finally, some of the studies that have 
argued against NFW haloes fitted their models in some ad 
hoc manner to the last measured data point or the maxi- 
mum of the rotation curve. However, there is no reason for 
treating a certain data point in a special way. In our analy- 
sis we fitted the models to all data points, weighted by their 
errors. 



6 SUMMARY 

We have analyzed high resolution HI rotation curves for 
a sample of 20 late-type dwarf galaxies. Taking beam- 
smearing and adiabatic contraction into account we have 
investigated to what extent these rotation curves put con- 
straints on the (central) density distributions of dark matter 
haloes. 

We have shown that two distinct degeneracies hamper 
a unique mass-decomposition. The first one, which we call 
the cusp-core degeneracy, owes to the fact that the observed 
rotation curves in general only sample the circular veloci- 
ties of the system at intermediate radii. No data is available 
at either very small or very large radii, making it virtu- 
ally impossible to discriminate between haloes with a con- 
stant density core and a r~ cusp. The other degeneracy, 
the mass-to-light ratio degeneracy, is well known from the 
rotation curves of (high surface brightness) spiral galaxies: 
as a result of the uncertainty in the stellar mass-to-light ra- 
tio, the relative amount of mass in the stellar disk can be 
exchanged with the amount of dark matter. It is noteworthy 
that several studies in the past have suggested that dwarf 
galaxies are not impeded by this mass-to-light ratio degen- 
eracy, and are therefore ideally suited to infer constraints 
on dark matter haloes. However, when the effects of beam 
smearing are taken into account, late-type dwarf galaxies 
are also plagued by the mass-to-light ratio degeneracy (see 
also the discussion in S99). 

In ~ 70 percent of the cases analyzed here, we find that 
the rotation curves are consistent with a ACDM cosmology. 
In the remaining ~ 30 percent, no meaningful fit to the 
observed rotation curves could be obtained with our mass 
models for any value of the inner slope of the halo density 
profile. Thus, although these galaxies are inconsistent with 
the ACDM model, they can not be considered to support any 
of the alternative dark matter models (unless a mechanism 



can be devised that can explain extremely small baryonic 
mass fractions) . This is most likely due to systematic errors 
and/or the fact that some of the assumptions underlying 
the models are incorrect. This emphasizes that care is to 
be taken when interpreting rotation curve fits; sometimes 
inconsistencies with CDM predictions are claimed without 
exploring the full freedom in halo parameters, and without 
addressing whether or not alternative models (i.e., with a 
constant density core) can yield realistic fits to the data. 

Our main conclusion, therefore, is that there is no con- 
vincing evidence against dark matter haloes of dwarf galax- 
ies having r _1 cusps. The HI rotation curves analyzed here 
are consistent with dark matter haloes with a = 1 and with 
concentrations as predicted for the currently popular ACDM 
cosmology. Together with the results for LSB galaxies pre- 
sented in BRDB and S99, we thus conclude that, based on 
the rotation curves of galaxies, there is currently no need to 
abandon the idea that dark matter is cold and collisionless. 
However, if future high resolution simulations confirm earlier 
findings of cusp slopes in the range of a = 1.5, or if it turns 
out that dwarf galaxies have disks with Tr 3> 1.0(M/L)q, it 
may be necessary to abandon CDM in favor of an alternative 
that yields haloes that are less steeply cusped. 

It is important to point out that the rotation curves 
studied here are also consistent with the presence of dark 
matter haloes with constant density cores. Thus, although 
current data does not require abandoning CDM, neither 
does it allow us to rule against its alternatives such as WDM, 
SFDM, or SIDM. Discriminating between these various dark 
matter models requires rotation curves of extremely high ac- 
curacy. Given the numerous sources for (systematic) errors, 
and the typical beam size of radio observations, we conclude 
that based on HI rotation curves alone at best weak limits 
on cosmological parameters and/or the nature of the dark 
matter can be obtained. In order to place more stringent 
constraints on the actual density distribution of dark mat- 
ter haloes one needs either data of higher spatial resolution, 
such as obtainable with Ha spectroscopy (i.e., Courteau 
1997; Blais-Ouellette, Carignan & Amram 1999; Swaters et 
al. 2000; Blais-Ouelette et al. 2000) or CO observations (e.g., 
Sofue et al. 1999) , and one needs to obtain independent con- 
straints on the disk's mass-to-light ratio from for instance 
stellar velocity dispersions (Bottema 1993) or, in the case of 
barred disk galaxies, from a detailed modelling of the veloc- 
ity fields (e.g., Englmaier & Gerhard 1999; Weiner, Sellwood 
& Williams 2001). 
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APPENDIX A: 
GALAXIES 



COMMENTS ON INDIVIDUAL 



UGC 731: As fo r th e models, which were moulded after 
this galaxy (see § 3.5), there is a large degeneracy in the 



model parameters. The only robust results seem to be that 
a 1.6 and /bar 0.05. Most importantly, the observed 
rotation velocities of this galaxy are in excellent agreement 
with CDM haloes, i.e., for a = 1, we find c ~ 16 — 2.3T_r. 
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UGC 3371: The best-fit models prefer a dark mat- 
ter halo with a steep central cusp. However, for a given Tr 
the minimum \ 2 ls achieved for a = a C rit, and these best- 
fit models are therefore unrealistic. Although \ 2 seems to 
depend rather strongly on a, this owes mainly to the ex- 
tremely small errorbars on V ro t- the best-fit rotation curves 
with T_r = 1.0 (M/L)0 and a = 0, 0.5, and 1.0 can hardly 
be discerned by eye. For a = lwe find c ~ 9.5— 1.45Th, and 
we thus conclude that UGC 3371 is consistent with CDM. 

UGC 4325: The quality of the fit improves consid- 
erably with increasing mass-to- light ratio up to T_r ~ 
7.5 (M/L)©, after which \ 2 increases rapidly. However, for 
Tr <; 7 (M/L)© the best-fit models have c = 1 and are 
thus unrealistic. Furthermore, S99 has shown that 0.5 ^ 
Tr 2.0 based on the stellar velocity dispersions, and 
Tr Jj 3 (M/L)q is unlikely in the light of stellar population 
models and UGC 4325's color of B - 7? = 0.85 (see Fig- 
ure [|). Clearly, the best-fit model is not the most realistic 
model. As is evident from the lower-left panel, models with 
Tr = 1.0 (M/L)© provide reasonable fits to the data, virtu- 
ally independent of a. For a = 1 we find c = 30.9 — 4.8T«, 
and we thus conclude that UGC 4325 is consistent with 
CDM. 

UGC 4499: The quality of the fit depends strongly 
on the stellar mass-to-light ratio. For Tr <; 2.0 (M/L)©, 
the models become unrealistic. The best fitting models have 
Tr = 0.0, which is also unphysical. For T_r = 1 (M/L)©, 
the best-fit models with a = 0.0, 0.5, and 1.0 yield virtually 
equally good fits to the data, but models with a ^ 1.1 are 
excluded by the data, since they require c < 1. For a = 1 
we find c ~ 9.0 - 7.0T fl . Henceforth, for Tr = 1.0 (M/L)© 
UGC 4499 is inconsistent with ACDM since it predicts a 
too small c. Consistency with B99's ACDM model requires 
Tr 0.7 (M/L)©. Since this is not unrealistic, we still con- 
sider UGC 4499 to be consistent with CDM. 

UGC 5414: The observed rotation velocities imply 
that Tr ^1.0 (M/L)©. Furthermore, the data favors a con- 
stant density core, and is clearly inconsistent with CDM. 
However, none of the mass-models provides a realistic fit, 
even for a = (either c < 1 or /bar < 0.01). 

UGC 6446: For the models to be realistic requires 
T R <, 2.5 (M/L)©. For T R = 1.0 (M/L)© the best-fit model 
has a ~ 1.5. However, this model is unrealistic (i.e., c = 1), 
but models with < a < 1.5 all provide virtually equally 
good fits. For a = 1.0 we find c ~ 17.4 - 10.3T fl and 
UGC 6446 is thus consistent with CDM. 

UGC 7232: The properties of this galaxy closely re- 
semble those of UGC 5414 and UGC 7323. Even for T R = 
do we not find a best-fit model with c > 1 and /bar > 0.01. 
We therefore conclude that no meaningful fit can be ob- 
tained for this galaxy. Note that the observed rotation curve 
consists of only five data points. 

UGC 7323: As for UGC 5414 and UGC 7232, no mean- 
ingful fit can be obtained for this galaxy. 

UGC 7399: In order for the models to be realistic re- 
quires Tr <, 4 (M/L)©. For a = 1 we find c ~ 23.1 - 7.0T R 
and UGC 7339 is thus consistent with CDM. 

UGC 7524: This galaxy has the best resolved rota- 
tion curve of all galaxies analyzed here. Unfortunately, equa- 
tion (^) does not yield a reasonable fit to the observed 
HI surface density, and we therefore opted to use the full- 
resolution data (see § ^ to model the unsmeared HI surface 



density distribution. The data favors T_r 4 (M/L)© and 
low values for a. However, for a = 1 we find c ~ 8.5 — 3.6T_r, 
consistent with CDM. 

UGC 7559: As is evident from the fact that Xrcd > 
18, none of our mass-models is able to yield a reasonable 
fit to the observed rotation curve. However, this is not too 
surprising since for r J; 0.9 kpc the velocity field is highly 
asymmetric (see S99) , which is not properly reflected by the 
errorbars. As for UGC 5414, UGC 7232 and UGC 7323 no 
meaningful fit can be obtained. 

UGC 7577: No results are plotted for this galaxy, since 
no model-fit was found to converge. We thus classify this 
galaxy as UGC 5414, UGC 7232, UGC 7323, and UGC 7603, 
in that no meaningful fit can be obtained. S99 suggested that 
UGC 7577, which is at a projected distance of only 37 kpc to 
NGC 4449, may be a dwarf galaxy that was formed by tidal 
interactions in the HI streamers around NGC 4449 (Hunter 
et al. 1998). Such tidal dwarf galaxies are expected to have 
little or no dark matter (Barnes & Hernquist 1992). 

UGC 7603: Similar to UGC 7524, no reasonable fit 
to the observed HI surface density can be obtained with 
equation (|l|) and we use the full-resolution data to model 
the unsmeared HI surface density distribution. As for the 
mass models: no meaningful fit can be obtained. 

UGC 8490: The rotation curve is well-resolved and 
prefers models with a low mass-to-light ratio. For a — 1 we 
find c ~ 24 — 10T_r, and this galaxy is thus consistent with 
CDM. 

UGC 9211: Models with T R = 1.0 (M/L)© provide 
reasonably good fits to the data, virtually independent of a. 
For a = 1.0 we find c ~ 18.3 — 3.6T_r, and we thus conclude 
that this galaxy is consistent with CDM. 

UGC 11707: This galaxy reveals a very large amount 
of freedom in its model parameters. This is partly due to 
the relatively large errorbars for the inner data points, and 
which is a reflection of the asymmetry between the receding 
and approaching rotation velocities at r Ss 7 kpc (see S99). 
For a = 1.0 one obtains c ~ 13.0 — 1.8T_r, and we thus 
conclude that this galaxy is consistent with CDM. 

UGC 11861: This is one of the few galaxies for which 
xtcd decreases with increasing Tr (see also UGC 5424, 
UGC 8490, and, to a lesser extent, UGC 7524 and 
UGC 12732). For Tr £ 3 (M/L)© the implied baryon 
fraction becomes unrealistically large. For a = 1 we find 
c ~ 16 — 3T_r, and we thus conclude that this galaxy is 
consistent with CDM. 

UGC 12060: The rotation curve of this galaxy is well 
fitted by the mass models. As for UGC 11707, there is a large 
amount of freedom in the model parameters. For a = 1 we 
obtain c ~ 34.6 — 9.8T_r. This implies that according to 
our definition UGC 12060 is inconsistent with ACDM for 
Tr = 0, since the best fit halo concentration is too large (i.e., 
c = 34.6 > c max = 31.9). However, for more realistic mass- 
to-light ratios the best-fit halo concentration is in excellent 
agreeement with the predictions of the B99 model, as we 
thus conclude that this galaxy is consistent with CDM. 

UGC 12632: The properties of UGC 12632 are similar 
to those of UGC 731. For a = 1.0 we find c ~ 14.0 - 1.3T«, 
and we thus conclude that this galaxy is consistent with 
CDM. 

UGC 12732: For this galaxy, larger mass-to-light ra- 
tios imply smaller values for a. For Tr JSJ 3.5 (M/L)© the 



© 2000 RAS, MNRAS 000, 



14 van den Bosch & Swaters 

resulting baryon fraction becomes unrealistically small. For 
T_r = 1.0 (M/L)q, the best-fitting model has a = 1.5, but 
also c = 1, and is thus unrealistic. For a — 1.0 we find 
c ~ 10.3 — 3.3T_r, and we thus conclude that this galaxy is 
consistent with CDM. 
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UGC 731 




0.5 1 1.5 2 2 4 6 8 

a radius [kpc] 

Figure Al. The panel in the upper-right corner plots the observed HI surface density (open circles with errorbars) together with the 
best-fit model of equation (|l|) both before (dashed lines) and after (solid lines) beam smearing. In the cases of UGC 7524 and UGC 7603, 
were no acceptable model fit can be obtained, the dashed lines indicate the observed HI surface density at the full resolution of the 
observations (see § which we use as a model for the unsmeared surface density of the HI. The four panels on the left show, from top 
to bottom, Xr e< j> c > ^200, and /t, ar of the best-fit models, as functions of a and for three different mass-to-light ratios (as indicated in the 
second panel). The lower-right corner panel, finally, plots the observed rotation curve (open circles with errorbars) together with four 
best-fit models with Yr = 1.0 (M/L)q: a = (solid lines), a = 0.5 (dotted lines), a = 1.0 (short-dash lines) and a = 1.5 (long-dash 
lines). These four models only differ in their dark matter properties; they have the same gaseous and stellar disks, the contributions of 
which are also indicated by (dot - short-dash) and (dot - long-dash) lines, respectively. These plots are useful for assessing the typical 
quality of the model fits and the cusp-core degeneracy discussed in the text. 
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Figure A3. Same as Figure Al but for UGC 4325 
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Figure A5. Same as Figure Al but for UGC 5414 
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Figure A7. Same as Figure Al but for UGC 7232 
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Figure A9. Same as Figure Al but for UGC 7399 
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Figure All. Same as Figure Al but for UGC 7559 
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Figure A12. Same as Figure Al but for UGC 7603 
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Figure A13. Same as Figure Al but for UGC 8490 
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Figure A15. Same as Figure Al but for UGC 11707 
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Figure A17. Same as Figure Al but for UGC 12060 
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Figure A19. Same as Figure Al but for UGC 12732 
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